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Abstract 

We demonstrate that sea quark effects of a magnitude expected from 
renormalization group considerations are clearly visible in the strong cou- 
pling constant measured in current full QCD simulations. Building on this 

(5) 
result an estimate of ajj^{Mz) is made employing the charmonium IS* — IP 

mass splitting calculated on full QCD configurations generated with two fla- 
vors of dynamical Kogut-Susskind quarks to fix the scale. 
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A distinctive feature of QCD that differentiates it from pfienomenological quark 
models of hadrons is the existence of sea quarks built into the theory. Nonetheless, 
finding a physical manifestation of sea quark effects has been an elusive subject in 
full lattice QCD simulations. In hadron mass spectrum calculations, for example, 
full QCD results for flavor non-singlet hadron masses agree with those of quenched 
QCD within statistical errors of 5-10% if the bare coupling constant for the latter 
is shifted by an appropriate amount [|T|. A similar situation holds for the critical 
coupling of the chiral transition at finite temperatures in full QCD; its value, 
though largely dependent on the sea quark mass, is reproduced quite well from 
that of the pure gauge theory by correcting for quark one- loop vacuum polarization 
effects §. 

A possible interpretation of the matching of full and quenched QCD by a shift 
of the bare coupling is that the shift represents an adjustment of the renormalized 
coupling constant at the low energy scale that dominates the behavior of quantities 
being simulatedp, ^. If this interpretation is valid, one expects that sea quark 
effects will become manifest in the renormalized coupling constant estimated for a 
scale sufficiently large compared to the dominant scale, since the full QCD coupling 
constant decreases more slowly than that of the pure gauge theory. In this article 
we present evidence that this in fact is the case: we find that the full QCD coupling 
constant extracted from two-flavor full QCD simulations is consistently larger than 
that of quenched QCD at large momenta ranging over fi ^ 3 — 7GeV when the 
scale is determined from the p meson mass; the difference in the magnitude of the 
full and quenched coupling constants is consistent with the picture that the two 
couplings merge when evolved down to the low energy scale /i < IGeV via the 
two-loop renormalization group. We also estimate the physical strong coupling 



constant for five flavors cnj^ at the Mz scale following the work of Ref, 



(5^ 
MS 

employing the IS* — IP mass splitting of charmonium states estimated for two- 



flavor full QCD. 

Calculation of the renormalized value of the coupling for a bare value ao taken 
in a simulation is facilitated by the recent study0 which has shown that lattice 
perturbative series is well convergent after lattice gluon tadpole effects are properly 
taken into account. A proposal for including tadpole effects in the relation between 
the bare and renormalized coupling is given by[0 

a^-ffin/a)-' = Pa^' + c, + NjCf + 0(«^), (1) 

where P is the plaquette expectation value, Cg = 0.30928 is the gluon one-loop 
contribution]^, and the last term represents the contribution of Nf flavors of 
quarks with c/ = —0.08848 for the Kogut-Susskind quark action[0 and c/ = 
—0.03491 for the Wilson actionjQ. Alternatively one may use the ay coupling 
defined from the static qq potential, which can be estimated from P via0 

- logP = yaif^)(3.41/a) (l + (d, + Ar^rf/)aSf'^ + 0(«^)) (2) 

with dg = —1.1855 and df = —0.0703 for the Kogut-Susskind quark action and 
df = —0.0249 for the Wilson action. The relation between the two couplings are 
given by[|, I 



,W-)/ N-l (Nf), ._i V , ri/^(^/)/,A2x 



1 



«7gtV)"' = ^vi^^r' + c^+0{a\^'^■\^^r), c^ = ^ (93 - lOAT^) . (3) 
Equivalently the A parameters are related by 

with bo = (11 - 2A^//3)/(47r)2. 

The renormalized coupling constant in the MS scheme extracted from (|l]) for 
quenched and two-flavor full QCD are compared in Fig. [l| for Kogut-Susskind and 
Wilson quark actions as a function of scale fi = it /a determined from the p meson 
mass. In full QCD the plaquette data are extrapolated linearly in the sea quark 



mass to rriq = 0. In quenched QCD we made a ninth order polynomial fit in (3 of 



plaquette values published in the literature [|14[ in order to calculate the values at 
/3 where data are not available. The trend is apparent in Fig. |l] that the full QCD 
coupling constant is systematically larger than that of the pure gauge theory when 
compared at the same scale /i. 

The solid lines in Fig. |1| illustrates the two-loop renormalization group evolution 
of the coupling constant. Deviation of a-j^ (n/a) from the solid lines toward 
smaller values of cutoff is in part ascribed to scaling violation effects due to a 
finite lattice spacing and in part to uncertainties of Ola"^) in the relation (|T]). One 
can estimate the magnitude of the latter through a comparison of the coupling 
constant extracted from (|l]) and (|^). This analysis shows that the latter estimate 
yields values for tt-p4 (^/^t) larger by about 3-5% at /i ^ 7GeV, and by about 
5-10% at yU ^ 3GeV. This is taken as uncertainties of our analyses. 

In Fig. 1^ we compare the two-loop renormalization group evolution of the 
full and quenched coupling constants toward small momenta /i < 0.5 — IGeV. 
The upper and lower edges of the bands in this figure correspond to ctjj^ (7r/a) 
estimated from the relation (|1|) and that from (^ including scale errors in order 
to take into consideration the two-loop uncertainty. For full QCD we employ the 
data taken at the highest (3 for the starting value, and for quenched QCD the one 
carried out at a value of (3 with a nearby value of // = ir/a. We observe that the 
evolution of the two coupling constants overlaps below /i ^ 0.4GeV, which is the 
dominant scale relevant for the p meson that is employed for fixing the scale. 

The results described above are fully consistent with the view that matching 
full and quenched results means adjusting the coupling constant at the relevant 
low energy scale, and that the full QCD couplings estimated for larger momenta 
should exhibit a slower decrease than the pure gauge coupling with a rate dictated 
by the renormalization group (3 function. 



Let us note that these findings provide support for the procedure of Refs. 0, 
^ for estimating the physical strong couphng constant from values measured in 
simulations with an incomplete spectrum of sea quarks. Namely one initially 
evolves the measured value at the cutoff scale down to the low energy scale /zq 
typical of the simulated hadron system using the flavor number of the simulation. 
The physical coupling constant Ug at that scale is equated to the evolved value, 
and as for larger scale is calculated through renormalization group incorporating 
the full spectrum of quarks active at each scale. A necessary condition for applying 
this procedure is that the scale /io is not too small in order not to spoil the two- 
loop approximation to the (3 function that breaks down for small momenta. The 
authors of Ref. [Q proposed to use the 15 — IP charmonium mass splitting for 
which potential models suggest /io ~ 0.4 — 0.75GeV. The advantage is that heavy 
quark propagators are easy to calculate and that the IS* — IP mass splitting is 
empirically insensitive to the quark mass, rendering its flne tuning unnecessary. 

We have carried out an analysis along this line employing the full QCD conflg- 
urations on a 20^ for two flavors of dynamical Kogut-Susskind quarks at (3 = 5.7 



with niqa = 0.01 [|T^. For charmonium spectrum measurement we used the Wil- 
son quark action for valence quarks |T7|, employing Gaussian smeared sources 
J2x,y'^x^ii^yf{x)f{y) and local sinks i'x^2i^x where Ti = Tg = f, 7^ a^'' for 
J/'ip, r]c, hc{op. 1) with /(x) = exp(— |a;p/4), and also Fi = 7^, F2 = cr^^ with 
f{x) = sin(2'7rxYL)/s(x) for hc{op. 2). We analyzed 72 conflgurations at K = 0.130 
and 75 conflgurations at i^ = 0.135 with the lattice size periodically doubled in 
the temporal direction. 

In Table |l] we list our result for the charmonium spectrum and the correspond- 
ing scale n/a extracted from the experimental value of the IP — IS* mass splitting 
457.8(5)MeV|15[ as input. The splitting is almost independent of the hopping 
parameter K which controls the charm quark mass, though it slightly depends 



on the choice of the operator for h^. [op. 1 or op. 2). The sphtting yields a value 
7r/a K, 7GeV, which is consistent with Ti/a = 7.01(28) [J16l estimated from the p 
meson mass. 

Our results for the physical strong coupling constant and the A parameter 
obtained with the scale listed in Table |I| are summarized in Table 0. In Table |^(a) 
the starting value is aj^{n/a) = 0.142 estimated from (|1|), while in Table ^(b) we 
use cty (3.41/a) = 0.169 obtained from (j^). For both cases the actual evolution is 
made in terms of the ay coupling since it is directly related to the heavy quark 
potential relevant for charmonium. We use /xq = 0.4 GeV and 0.75 GeV for the 
matching scale and take an average of the two choices for the central value for the 
strong coupling constant. The errors are estimated by allowing /xq to vary over the 
range /xq = 0.4 — 0.75GeV and the scale n/a within the quoted error. 

Our results are consistent with the previous lattice estimates carried out in 
quenched QCD^, Q. Compared to the world average of phenomenological deter- 
minations ajj^{Mz) = 0.118(7) [|1^, the results are somewhat small especially for 
those estimated from the relation ([1|) . We do not view the difference to be alarm- 
ing at this stage since there exists a 5% uncertainty in our value of ajj^{Mz) due 
to that of the input value at the cutoff scale, and an additional 5% that results 
from the matching procedure, as well as errors in the experimental value. 

To summarize, our analyses have shown that sea quark effects of a magnitude 
expected from renormalization group considerations are visible in the strong cou- 
pling constant measured in current full QCD simulations incorporating up and 
down quarks. This indicates a promising prospect for a realistic determination of 
the strong coupling constant including the full spectrum of sea quarks since incor- 
porating heavy quarks such as strange and charm is not difficult from the view of 
the necessary computer power compared to that for light quarks. 
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K rriri^a 



m 



J/i' 



,a 



mua 



Arriip^isa 7r/a(GeV) 



0.130 1.474(16) 1.466(17) op. 1 

op. 2 

0.135 1.285(16) 1.275(15) op. 1 

op. 2 



1.692(53) 0.221(49) 6.5(1.5) 

1.672(32) 0.200(32) 7.2(1.2) 

1.504(53) 0.221(48) 6.5(1.5) 

1.487(32) 0.205(35) 7.0(1.2) 



Table 1: Charmonium spectrum for Wilson valence quarks in full QCD at /3 = 5.7 
with two flavors of dynamical Kogut-Susskind quarks with niqa = 0.01. 



K 


a^(5GeV) 




(y§s{rnz) 


^(5) 
MS 


(a) estimate from aii^(7r/a) 






0.130(oj9. 1) 


0.168t°:°l^ 


1411II MeV 


n 1 04+0-005 


921^1 MeV 


0.135(op. 1) 


0.168t°:°lg 


14lli 


n 1 04+0-005 


92t^2 MeV 


0.130(op. 2) 


0.173i°:°^^ 


159l^t 


0.1061°:°°^ 


105t^^ MeV 


0.135(op. 2) 


n 1 vo+0-012 

U-J- 1^-0.009 


1541^? 


o.io6t°:°°^ 


102ti MeV 


(b) estimate from a^ (3.41 


/a) 






0.130(op. 1) 


0.1861°:°}^ 


2QQ+_ll MeV 


r. 1 1 Q+0.006 
U-J-iU_o.005 


142t^^ MeV 


0.135(op. 1) 


o.i86J:°:°li 


20611^ 


n 1 1 n+0.006 
U--LiU_o.005 


142ti MeV 


0.130(op. 2) 


0.1921°:°!? 


23lttl 


n 1 1 O+0.004 
U.iiZ_Q004 


16lt^^ MeV 


0.135(op. 2) 


0.1901°:°}^ 


224t^° 


n 1-19+0.004 
U.ilZ_oo04 


1561^ MeV 













Table 2: Strong coupling constant and A parameter in the MS scheme according 
to the Particle Data Group definition calculated with the IS* — IP charmonium 
mass splitting with Wilson valence quarks for fixing the scale. 



Figure captions 



Figure 1: Comparison of aij^ (vr/a) estimated via (P for two-flavor full QCD 
(filled symbols) and quenched QCD (open symbols) with the scale fixed by the p 
meson mass. 
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Figure 2: Evolution of a-p-i- (/i) for quenched and two- flavor QCD. Bands cor- 
respond to uncertainties in the estimate of cnjj^ {jj). Arrows indicate the start- 
ing value taken from (a) Fukugita et al. {Nj = 2, /? = 5.7)|jl0[ and Sharpe et 
al.{Nf = 0,/? = 6.2) ini, and (b) Gupta et al.{Nf = 2, P = 5.6)0 and Butler et 
aliNf = 0, P = 6.17)^. 
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This figure "figl-l.png" is available in "png" format from: 
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Fig. 1(a) 
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Fig. 1(b) 
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(b) 1/ oc— (|i) Wilson quark action 
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